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Abstract 

SU(2) lattice gauge theory is extended to a larger coupling space 
where the coupling parameter for horizontal (spacelike) plaquettes, 
Ph-, differs from that for vertical (Euclidean timelike) plaquettes, f3y. 
When /3jj — > oo the system, when in Coulomb Gauge, splits into 
multiple independent 3-d 0(4) Heisenberg models on spacelike hy- 
perlayers. Through consideration of the robustness of the Heisen- 
berg model phase transition to small perturbations, and illustrated 
by Monte Carlo simulations, it is shown that the ferromagnetic phase 
transition in this model persists for (3h < oo. Once it has entered the 
phase-plane it must continue to another edge due to its symmetry- 
breaking nature, and therefore must necessarily cross the j3y = fin 
line at a finite value. Indeed, a higher-order SU(2) phase transition 
is found at /3 = 3.18 rb 0.08, from a finite-size scaling analysis of the 
Coulomb gauge magnetization from Monte Carlo simulations, which 
also yields critical exponents. An important technical breakthrough is 
the use of open boundary conditions, which is shown to reduce system- 
atic and random errors of the overrelaxation gauge-fixing algorithm by 
a factor of several hundred. The string tension and specific heat are 
also shown to be consistent with finite-order scaling about this critical 
point using the same critical exponents. 



1 Introduction 



It would be very nice to have a local order parameter for confinement. On 
large lattices the Polyakov loop becomes extremely small in both phases, 
as do large Wilson loops. In addition, it is difficult to accurately split off 
the area-law confining term of Wilson loops from the other terms present. 
For spin models with spontaneously broken symmetries it is much easier to 
study the ferromagnetic phase transition, due to the local magnetic order 
parameter which (along with its moments) is easy to measure on any size 
lattice. 

Indeed, there actually is such a local order parameter, the ferromagnetic 
transition of which signals deconfinement. It requires, however, fixing con- 
figurations to minimal Coulomb gauge. Since this is an iterative gauge-fixing 
procedure, it is fairly costly, but the benefit of a local order parameter may 
be worth it. Nearly two decades ago, it was envisioned that this order param- 
eter could become a prime means for studying confinement [T] . Difficulties in 
gauge-fixing have largely prevented these techniques from becoming main- 
stream. 

The minimal Coulomb gauge involves maximizing the trace of all links 
in three of the four lattice directions (these three are called "horizontal" di- 
rections below). The fourth-direction (vertical) links can be reinterpreted as 
0(4) "spins," the local order parameter. The minimal Coulomb gauge has a 
set of remnant SU(2) symmetries, global on the three horizontal-directions, 
but still local along the fourth. This is because if links are written as ao+ia-T 
the do components of the horizontal links are invariant under such a gauge 
transformation, and this is the quantity being maximized to satisfy the gauge 
condition. Separate SU(2) transformations on the two ends of the fourth- 
direction links, which are being treated as spins, generate for them a global 
SU(2) xSU(2) = 0(4) symmetry, which breaks spontaneously if the spins 
magnetize. In the limit f} — > <x>, the gauge condition sets the horizontal links 
to unity, and the plaquette interaction collapses into a spin interaction, a 
dot-product between the 0(4) spins. Thus, in this limit, the SU(2) lattice 
gauge theory becomes a set of 3-d 0(4) Heisenberg models (or non-linear 
sigma models) at zero temperature (using T = as temperature). As 
/3 becomes finite, then the primary interaction is still the Heisenberg one, 
since the horizontal links still are on average close to the identity, but the 
deviations from the identity introduce new interactions between the spins, 
which vary spatially. For large (3 these interactions are either small or, if 
large, rare. There are no direct interactions between vertical links in dif- 
ferent hyperlayers; they interact indirectly through shared horizontal links. 
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For a fixed set of horizontal links, the spin-model now has the appearance of 
a spin-glass, with the deviations of horizontal links from the identity intro- 
ducing spatially-varying "disorder interactions." As /3 is lowered, both the 
pseodo-temperature, T — 1//3, and the amount of disorder from the deviation 
of horizontal links from the identity increases. Thus, one expects that the fer- 
romagnetic phase that must exist at = oo (since the 3-d 0(4) spin model is 
ferromagnetic at zero temperature), should either have a transition to a para- 
magnetic phase, or possibly first to a spin-glass ordered non-ferromagnetic 
phase followed by a second transition to a paramagnetic phase. Because the 
strong coupling limit is completely random, the system must eventually enter 
a paramagnetic phase. 

Conventional thinking would have the transition to paramagnetic phase 
happen immediately at T = 0+, putting the entire Wilson axis except /3 = 
oo in the paramagnetic phase. However, ferromagnetism in the Heisenberg 
model is quite robust against the introduction of disordering interactions. 
For instance the ±J 3-d 0(3) spin-glass model remains ferromagnetic at 
zero temperature until 21% of randomly selected interactions are changed 
from ferromagnetic to anti-ferromagnetic [2J. The 0(4) model is expected to 
behave similarly. Thus it is a little hard to imagine how the infinitesimally 
small disordering interactions introduced at T = 0+ in the gauge theory 
can so effectively kill the ferromagnetic order. Recently, using spin-glass 
methods, a spin-glass to paramagnet transition was observed at = 1.96, 
strongly supporting the second hypothesis above[3]. Unlike the ferromagnetic 
phase, which is necessarily unconfined [I], the spin-glass phase, which has a 
hidden pattern of frozen disorder, is still confining. Below, the necessary spin 
glass to ferromagnet phase transition is searched for. That itself could be at 
T = 0+, again consistent with the whole /3-range being confining, however, 
the evidence presented below from a variety of angles points to this transition 
taking place on the infinite lattice around /3 = 3.2. 

A previous attempt at finding this transition showed a possible transition 
around /3 = 2.6[5J. However this was not consistent with the apparently 
obvious confinement in the region /3 = 2.7 to 2.85 from the Polyakov loop 
and interquark potential^ U\- It was later discovered that the gauge fix- 
ing procedure was not working well during Polyakov-loop tunneling events. 
This depressed the magnetization on smaller lattices resulting in a prema- 
ture crossing of the Binder cumulant. It is not surprising that the local 
gauge-fixing algorithm has trouble with the global constraints imposed by 
the gauge-invariant Polyakov loops. Switching to open boundary conditions 
seems to completely eliminate this pathology, and as a result not only lessens 
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the systematic error, but vastly reduces the random error in quantities such 
as the Binder cumulant (by a factor of several hundred). Evidence is given 
below that the remaining systematic error is at most of order, and probably 
much less than, this much-reduced random error, which itself is 100 times 
smaller than the size of ordinary thermal fluctuations between configura- 
tions. Even with the reduction of random error that results from averaging 
over gauge configurations, the possible systematic error remains less than or 
of order the net random error. Since the crossings identifying the new phase 
transition are verified at several /3's to more than 9a, systematic error no 
longer appears to be an issue. Indeed, the use of open boundary conditions 
seems to have the potential to transform Coulomb-gauge studies to the high- 
precision realm. The use of open boundary conditions in gauge theories has 
recently been justified by Liischer and Shaefer[8] and used successfully by 
them to reduce barriers to changes in the topological charge. The negative 
features of open boundary conditions, such as the lack of translational invari- 
ance, are not so important for the large lattices in use today, and are rather 
easily dealt with. 

In section 2 the effectiveness of the Coulomb Gauge relaxation using open 
boundary conditions is detailed. In section 3, the larger coupling space where 
purely horizontal plaquettes have a different coupling parameter, (3h, than 
those with a vertical link, /3y, is explored. This plane includes both the 3-d 
0(4) Heisenberg model on the 0h — 00 edge, and the ordinary 4-d SU(2) 
gauge theory along the (3h — 0v line. A transition very similar to the 
Heisenberg transition in terms of critical exponents is seen for large but non- 
infinite Ph- At this position the transition has not so high critical exponents 
as seen for the fiv = Ph case and is very easily studied with standard methods 
(Binder cumulant crossings etc.). An important point is that once it has 
been established that the transition enters the ((3h,Pv) coupling plane from 
the Heisenberg edge, it must continue until it hits another edge. Due to 
the symmetry-breaking nature of the transition, the coupling plane must 
be split into separated symmetry-broken and unbroken regions. As (3h is 
reduced, the critical f3y grows, thus a transition at finite (3 = (3y = Ph has 
to exist. Because the large (3h transition is so similar to the Heisenberg 
transition, and so easily established with standard Monte Carlo methods, 
this argument is perhaps the most solid evidence yet that zero-temperature 
SU(2) lattice gauge theory must have a phase transition at finite /3. This 
transition breaks the remnant symmetry remaining after the Coulomb gauge 
fixing. Such a ferromagnetic phase is known to be non-confining [3]. In short, 
from a statistical mechanics point of view, what is being shown here is that 
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SU(2) lattice gauge theory, like many other lattice gauge theories, has a lot 
in common with spin theories with the same symmetry once the gauge is 
fixed. It is essentially a magnet. 

The SU(2) transition itself is studied with Monte-Carlo simulations in 
section 4. The open boundary conditions allow for much more precise 
Coulomb-gauge data than previously available. Here the transition, seen 
at ft — 3.18 ±0.08, is somewhat harder to study from a technical standpoint, 
because the critical exponents and the cumulant crossing values are fairly 
extreme (in comparison to the Heisenberg model). However, even here the 
cumulant crossings can be verified to more than lOer. 

Because the effective 0(4) spin models are on 3-d hyperlayers, there are 
some interesting dimensional issues concerning critical behavior. When these 
are sorted out, fits give a correlation- length exponent v — 1.7 for the 4-d 
gauge theory. This makes the corresponding singularity in the specific heat 
much too soft to be seen numerically. However, the specific heat (which 
shows almost no finite lattice size dependence) can still be fit to the infinite- 
lattice functional form associated with such a transition. In section 5 it is 
found that the range of v to which the sharp rise of specific heat in the 
crossover region can be fit is consistent with this value, as is the v found 
from a similar fit to string tension measurements and from /3-shift fits to 
Polyakov-loop transition points for asymmetric lattices. Thus the behavior 
of these important gauge-invariant quantities is consistent with correlation 
length scaling associated with a higher order phase transition at ft — 3.2, 
using the critical exponent found from the Coulomb-gauge magnetization 
scaling behavior. The final picture that emerges is quite unconventional. The 
SU(2) lattice gauge theory appears to have three phases: a strong-coupling 
paramagnetic phase from ft = to 1.96, a spin-glass phase from 1.96 < 
ft < 3.2 and a non-confining ferromagnetic phase for ft > 3.2. The existence 
of the spin-glass phase explains why there is a gap (the crossover region) 
between regions that are well-approximated by strong and weak coupling 
approximations. 

2 Coulomb gauge fixing and boundary con- 
ditions 

The first thing to be mentioned here is that these are not technically 
"Coulomb gauge simulations." To actually stay within the Coulomb gauge 
during Monte Carlo simulation would be difficult and would also require 
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adding ghost fields to take care of the Fadeev- Popov determinant. Rather, 
the Coulomb gauge fixing here is part of the post-simulation measurement 
procedure, applied to a completely conventional gauge-invariant Monte-Carlo 
simulation. In fact, provided that the gauge fixing algorithm works well, the 
absolute- value of the resulting fourth-direction link magnetization is actually 
gauge-invariant in the following sense. If a gauge configuration is subjected to 
a random gauge transformation and then reprocessed through the Coulomb 
gauge fixing algorithm, the same result should ensue, provided there is a 
unique minimum. The direction of the magnetization will be different due 
to the possible inclusion of a remnant global symmetry transformation, but 
all important results are dependent only on the absolute value of the mag- 
netization and its even moments, as usual. The minimal Coulomb gauge 
on a finite lattice is not generally subject to the Gribov ambiguity, because 
the chances of not having an absolute minimum are quite small. However, 
it is also numerically difficult to find an absolute minimum and relaxation 
algorithms will likely find local minima. So, whereas in principle there is a 
unique solution, in practice many solutions are found and in that practical 
sense there is still a "Gribov problem" on the lattice. 

The method used here to set the Coulomb Gauge is iterative overrelax- 
ation, with a 70% overshoot. There are a number of techniques which have 
been tried to improve on this basic algorithm. One is to subject the initial 
configuration to random gauge transformations and then relax each of these 
separately. This is done, say, 10 times. The one which achieves the best 
maximization of the traces of horizontal links is then chosen, and the others 
discarded. Another more-sophisticated approach is simulated annealing[9|. 
Both are quite demanding of computer time. 

The basic overrelaxation algorithm utilizes local gauge transformations 
at each site, cycling through the lattice multiple times. From 40-4000 sweeps 
through the lattice are needed for reasonable convergence (the larger number 
for lower (3). However, with periodic boundary conditions, this does not 
exhaust all of the symmetry transformations available. There are still the 
Z(2) 4 Polyakov loop transformations in which a hyperlayer of spins pointing 
in a particular direction are all multiplied by -1, the non-trivial element 
of the center. It has been found that adding these transformations in a 
preconditioning step where the best of the 16 Polyakov-loop sectors is found 
first, before overrelaxation, improves the final results considerably |10] . 

When one uses multiple relaxations, most give similar results, however 
now and then a really bad solution is found (local minimum far from the 
global one). So it is likely that using only a single relaxation would result 
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in too many outliers; using more relaxations will reduce these but at an 
obvious computational cost. The technique chosen for the current study is 
somewhat different. The idea is to couple the gauge- fixing more closely with 
the original simulation. Gauge fixing is performed after each sweep of the 
simulation, which uses a Metropolis algorithm with about 50% acceptance. 
In most cases, new gauge elements are not chosen from the full group, but 
from a limited neighborhood of the original link. Throughout most of this 
study, this neighborhood was chosen to be the nearest half of the group space, 
but the adjustability of this hit-size also gives one a handle to test the effi- 
cacy of the gauge-fixing algorithm. The previous gauge-fixed configuration 
is used as the starting point for the Monte Carlo sweep. Thus, since only 
half of the links are changed, and those that do are mostly changed by small 
amounts in each sweep, the gauge is still "partially fixed" at the beginning 
of the next overrelaxation. This better starting point seems to avoid the 
"really bad minima" which are sometimes found when starting from a ran- 
dom gauge transformation. Since the overrelaxation is done even during the 
initial equilibration portion of the simulation, it has a long time to find a 
good "groove" and maintain it while the gauge configuration slowly evolves. 
Although the preconditioning step mentioned above is included as well, it 
turns out not to work with this close-coupled algorithm, because, since af- 
ter only one sweep the memory of the previous relaxation is still there, the 
algorithm always prefers not to jump to another Polyakov loop sector. One 
could possibly spawn 16 new configurations and relax each one, but that is 
simply too costly. 

This algorithm seems to perform well most of the time, but it was no- 
ticed that at high-beta (f3 > 2.9), where the magnetization was quite large 
and histograms peaked very far from zero, occasional "storms" of low values 
occurred. This produced a small secondary peak about zero magnetization 
(Fig. 1-lower). At first this was possibly thought to be a real phenomenon, 
and evidence for a first-order phase transition, but the persistence of this 
second peak over a large range of f3 was not consistent with that interpreta- 
tion. It was later discovered that the observed "storms" of low magnetization 
appeared to occur during Polyakov loop tunneling events. The known sensi- 
tivity of the gauge-fixing algorithm to Polyakov loop sector changes pointed 
to it as the cause of the storm - it was simply having trouble finding a good 
minimum during the tunneling. Since this represents a global shift in the 
vacuum, akin to a quantum phase transition, it is not surprising that the 
local gauge- fixing algorithm has trouble initially following it. After a few 
tens of sweeps it finds a new groove, but even a few percent of bad sweeps 
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Figure 1: Coulomb gauge magnetization histograms at j3 — 3.5 on 16 4 lat- 
tices for periodic boundary conditions (lower) and open boundary conditions 
(upper). The lower peak on the lower graph contains very few actual config- 
urations, bars being magnified by the l/|m| 3 geometry factor. Left scale is 
for lower graph and right for upper. 



with very different values will affect moments of the magnetization unaccept- 
ably. The memory effect of the close-coupled algorithm is a negative feature 
when it comes to following sudden vacuum tunnelings. These problems only 
appear to occur at values of /3 above the point at which the Polyakov loop 
is broken on a given lattice (e.g. (3 > 2.7 on the 16 4 lattice). 

In order to test this further, the algorithm was tried with open boundary 
conditions (OBC). Results were beyond all expectations. The magnetiza- 
tion distribution completely loses the suspicious second peak (Fig. 1 - upper 
graph), and the gauge condition function no longer has large excursions. 
These reduced fluctuations also resulted in much lower statistical errors on 
quantities such as susceptibility and Binder cumulant. This suggests that 
much of the observed error in the periodic boundary condition (PBC) simu- 
lations did indeed come from fluctuations in the efficacy of the gauge setting 
algorithm. 

The magnetization order parameter is defined first by defining an 0(4) 
vector a = (a , a±, a 2 , a 3 ) from each 4th-direction pointing SU(2) link a l + 
i Y^=i a j T j an d averaging these over each 3-d perpendicular hyperlayer. 

1 

m=jj ]T a. (1) 

hyperlayer 
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Figure 2: Results of multiple gauge relaxations on the same gauge config- 
uration after random gauge transformations for (a) periodic boundary con- 
ditions and (b) open boundary conditions. The single square datapoint on 
each graph is the result from the closely coupled algorithm. Trend lines in- 
dicate the remaining degree of correlation between magnetization and gauge 
function. 



Expectation values are then taken as ensemble averages over both hyperlay- 
ers and gauge configurations. In Fig. 2, scatter plots are shown for multiple 
relaxations following random gauge transformations on the same gauge con- 
figuration, one taken from a 16 4 lattice at /3 = 3.5. (Note it is the high f3 
region that is problematic in the PBC simulations and data in this region is 
also the most important for this paper). The gauge function being maximized 
is shown on the horizontal axis, and the magnetization |m| on the vertical 
axis with (a) showing the PBC case and (b) the OBC case. Note the large 
difference in scales. The scatter for the OBC case is seen to be smaller by 
a factor of about 200. Also for OBC there is very little correlation between 
the small remaining fluctuations in gauge function with the observable of 
interest, the magnetization. This would seem to indicate that the residual 
error is mostly random, rather than systematic, and thus easily observed. 
This contrasts with the rather strong correlation seen with PBC, which sug- 
gests a systematic error at least as large as the random. In addition, the 
OBC fluctuations are around 100 times smaller than the ordinary thermal 
fluctuations in magnetization between different gauge configurations, adding 
less than one part in 10 4 to the susceptibility, which is completely negligible 
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in comparison to the ordinary statistical errors in this study. Thus, switch- 
ing from PBC to OBC would seem to transform use of Coulomb gauge with 
Monte Carlo data from a relatively rough and problematic method into a 
precise one, where the additional spurious fluctuations in quantities due to 
imprecision in gauge fixing are negligible compared to the ordinary thermal 
(quantum) fluctuations. 

Another test used for finding systematic error was to vary the Metropolis 
hit size for the run, using the closely-coupled algorithm where the gauge 
relaxation is performed after each sweep. This is the maximum amount by 
which new gauge elements are allowed to differ from old ones in an update. 
Hit sizes vary from -1 to 1, with -1 meaning no change allowed at all and 1 
meaning the full group (it is the negative of the minimum a component of the 
update matrix). One would expect that the gauge setting algorithm would 
perform better when smaller changes were being made, since it would be less 
likely to find a local minimum far from the global one, and basically have more 
tries on a more restricted set of configurations. These runs were compared 
with those using a single relaxation after a random gauge transformation, 
and also the best of five and best of ten relaxations after random gauge 
transformations. The latter were only run at a hit size of zero because they 
would not be expected to show a dependence on hit size. The results are 
shown in Fig. 3 for Binder cumulant and magnetization. It is clear that 
there is very little difference between any of these algorithms - even the single 
relaxation after random gauge transformation gives results consistent with 
the others. For the magnetization, a slightly larger value is obtained with 
best of ten method than for a single relaxation. The best of ten value agrees 
with that from the closely-coupled algorithm. This suggests that the closely- 
coupled algorithm with zero hit size, as used in this paper is competitive with 
the best of ten algorithm. However the differences are so small it appears 
that even the single relaxation after random gauge transformation could be 
used in the OBC case. These runs were long enough (50,000 sweeps) that 
they should have encountered exceptional configurations if such still exist, 
one or two bad values from which would have had a large effect on the 
higher moments. So the data seem to indicate that the previous difficulties 
with Coulomb gauge relaxations no longer occur when using open boundary 
conditions. Indeed using open boundary conditions opens the path to highly 
precise measurements in the Coulomb gauge in which the efficacy of the gauge 
fixing algorithm is no longer a factor. 
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Figure 3: Binder cumulant, U, and magnetization vs. hit size for the closely- 
coupled algorithm on 16 4 lattices at (3 = 3.5 (diamonds). Also shown is 
relaxation after a single random gauge transformation (open circle), best of 
five (x) and best of 10 (filled circle). Points are horizontally offset slightly 
at hit size for clarity. Also shown are 24 4 values for U (triangles), demon- 
strating that algorithmic differences, if they exist, are much smaller than the 
gap between 16 4 and 24 4 values. 
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3 Expanding the parameter space to connect 
the spin model to the gauge theory 



The possible existence of a zero-temperature phase transition in the usual 4-d 
SU(2) gauge theory can be put into clearer focus by considering a larger cou- 
pling space in which horizontal and vertical plaquettes have different coupling 
parameters, (3h, fiv- Vertical plaquettes are those that include the fourth 
(Euclidean time) direction. As discussed above, for 0h = oo horizontal links 
are locked to unity in the Coulomb gauge, which transforms the 4-link gauge 
interaction into a 2-link nearest-neighbor spin interaction (scalar product) 
among the vertical links. Thus the theory is exactly a set of multiple non- 
interacting copies of the 0(4) 3-dimensional Heisenberg model (the quantum 
version of which is the non-linear sigma model). A similar connection of the 
SU(2) gauge theory to the Heisenberg spin model in one fewer dimensions was 
explored analytically in [TTJ. The 3-d 0(4) classical Heisenberg model has 
a well-known higher-order ferromagnetic phase transition at /3y = 0.936(1) 
with critical exponents v = 0.7479(90), 7/1/ = 1.9746(38), fi/v = 0.5129(11), 
and a = 2 — dv = —0.244(27) [12]. Since a < there is no infinite sin- 
gularity in the specific heat, but rather a cusp. As mentioned before, the 
long-range ferromagnetic order of the Heisenberg model is quite robust to 
the addition of a significant number of randomly chosen antiferromagnetic 
links, deleted links, or added next neighbor interactions j2j [13]. (Although 
the 0(3) model is usually studied, the 0(4) model appears to differ only in 
detail.) Thus, it would be quite surprising if the transition did not enter the 
WviPh) coupling plane as (3h is backed off from 00. In other words one would 
expect the transition to still exist for large but finite flu- In this region, the 
horizontal gauge fields will still be very close to unity, so the Hamiltonian 
would still be primarily the Heisenberg nearest-neighbor interaction. Small 
additional terms of order 1//3h would be present which would add disorder 
and likely add next-neighbor interactions if the horizontal gauge fields were 
to be integrated out. A renormalization of the pseudotemperature (l//3y) 
would also be expected. Because the integrated horizontal gauge fields do 
not break the remnant 0(4) symmetry which remains after Coulomb gauge 
fixing, these effective spin interactions cannot favor domain formation as a 
random external magnetic field would[14j. If the phase transition does enter 
the coupling plane for (3 H < 00, then due to its symmetry-breaking nature 
it cannot end except at another edge of the phase-plane (it must divide the 
phase-plane into two disconnected regions of different symmetry). As 0h is 
lowered from infinity, the critical value of /3y must increase to compensate 
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for the additional disorder, so clearly these must then meet at a finite value 
Pc — Pvc = 0h, which would be the critical point of the SU(2) gauge theory. 

Since a ferromagnetic phase in Coulomb gauge is necessarily 
deconfined[l], this calls into question the customary assumption of confine- 
ment in the continuum limit (f3 — > oo). If the conventional interpretation of 
no phase transition in SU(2) is correct, then the 0(4) critical point would 
have to exist only for (3 H = oo and cease to exist for any fi H < oo. In other 
words the slightest 1/(3h > coupling must destroy the long-range order. 
Since the large (3h case turns out to be simpler and more easily studied than 
the fly = Ph symmetric case, it seems worth focusing on it first to investi- 
gate this important question. It might, for instance, be possible to extend 
the proof of long range order in the Heisenberg model [TJ] to the small 1//3h 
region, which by the above argument would prove the existence of an SU(2) 
transition. 

A Monte Carlo simulation at (3h = 20 was performed to look for a 
Heisenberg-like transition in the large 0h region. Such a transition is eas- 
ily found using the Coulomb gauge method with open boundary conditions. 
As previously mentioned, the order parameter is simply the fourth-direction 
pointing links which act as 0(4) spins. The "gold standard" method of find- 
ing a Binder cumulant crossing from simulations on different size lattices is 
used to find the position of the phase transition on the infinite lattice, just 
as is done for the case of the spin model itself. The Binder cumulant [16J, 
U = 1— < |m| 4 > /(3 < |m| 2 > 2 ), is plotted in Fig. 4, showing a definite 
crossing near (3y = 1.01. 

For the 0(4) order parameter, Binder cumulant values range from 0.5 in 
the deep paramagnetic region to 2/3 in the deep ferromagnetic region[17J. 
Fig. 5 shows scaling plots for Binder cumulant, U, subtracted susceptibility 
X = V(< \fn\ 2 > — < \m\ > 2 ), and magnetization < \m\ >, demonstrating 
good collapse to finite size scaling ansatze [18] (see axis labels for detailed 
scaling functions, with T = 1/ f3y and T c = l//3y c ). 

One disadvantage of using open boundary conditions is that lack of 
translational invariance makes the Fourier transform impossible, and thus 
the second moment correlation length which also yields scaling in- 
formation, cannot be defined. However, the unsubtracted susceptibility, 
Xb = V < \m\ 2 > when raised to the appropriate power can be used as 
a surrogate for the correlation length, 

U = X b h - (2) 

The subscript L refers to the linear lattice size. The quantity C,l/L can 
be then used in the collapse fits to further constrain the determination of 
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Figure 4: Binder cumulant crossing for the (3h = 20 theory. U24 exceeds U 2 q 
by 2.5a and Uiq by 9cr at /3 = 1.06. Diamonds are 16 4 , squares 20 4 , and 
triangles 24 4 lattices on this and all subsequent plots. 

exponents 7 and v (see Fig. 5d). Once 7/V has been determined, Xl/L may 
be used to see crossings at a fixed point corresponding to the infinite lattice 
critical point. Fig. 6 gives a very clear crossing at a coupling consistent 
with that of the Binder cumulant. The overall fit giving the best collapse 
of the four quantities, shown in Fig. 5, gives (3v c = 1-01(2), v = 0.97(17), 
j/u = 2.00(25), 0/u = 0.51(14), with an overall x 2 /d-f- = 2-4. Only the 20 4 
and 24 4 data are used in the scaling fits as the 16 4 appears a bit too small to 
follow the universal scaling functions within the statistical errors here (higher 
order corrections are sometimes needed for smaller lattices). These runs were 
all for 50,000 sweeps following 10,000 equilibration sweeps. The 7/1/ and /3/u 
values are consistent with those for the 0(4) Heisenberg model, whereas the 
value of v appears somewhat larger. This makes sense considering that the 
value of v seen below for the 0h — Pv SU(2) case is much higher, around 3.42. 
Thus the /3# = 20 transition appears to have already started evolving in this 
direction but is otherwise very similar to the Heisenberg transition itself. Of 
course the main result here is that the 0h = 20 transition simply exists at 
a finite fiy c . As argued above, due to the symmetry breaking nature of the 
phase transition, this is a sufficient condition for the existence of an SU(2) 
transition at finite f3. Binder cumulant crossings and scaling collapse plots are 
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Figure 5: Finite size scaling collapse plots for U, x, < l^l >> an d £l/L 
the p H = 20 theory. 
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Figure 6: Crossing graph for for (3 H = 20. The 24 4 value exceeds the 
20 4 by 21a and the 16 4 by 43a at = 1.06. 

tried and true methods for finding phase transitions. Calling into question 
the (3h = 20 result would also seem to call into question the Heisenberg 
results themselves, since the methods and data are so similar. The reliability 
of these methods has been demonstrated in a large variety of models. If 
something drastically different were to happen on very large lattices, it is 
unlikely that such clean crossings would have been observed here. In order 
for the transition to disappear, the behavior of U and C,l/L in the weak 
coupling region would have to switch from the observed increasing function 
with lattice size to a decreasing function. Such a change in behavior would 
be highly unusual, if not unprecedented. 

The suggested phase diagram in the extended plane is given in Fig. 7. 
The Heisenberg transition is seen on the upper boundary with the 0h — 
20 transition, just described, below it. Further down, along the (3y = Ph 
symmetry line, is the SU(2) transition at (3 = 3.2 to be described in the 
next section. Also shown is the spin glass transition for SU(2) at = 1.96 
described in [3] . The point where the spin-glass phase first splits off has not 
yet been determined, nor has the behavior below the symmetry line. The 
bottom boundary is a set of one-dimensional 0(4) spin models, which are in 
the unbroken phase except at fiy — oo. So a likely place for the transitions 
to end up is this point, but further investigation is needed. 
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Figure 7: Suggested phase diagram from the phase transitions seen in this 
paper at /3h = 20 (square) and (3h = Pv SU(2) gauge theory (triangle). Also 
shown are the 3-d 0(4) Heisenberg model transition on the upper axis (star) 
and the paramagnet to spin-glass transition seen in ref.[3] (diamond). 

4 Monte Carlo simulations 

A similar study was performed for the regular SU(2) theory, again using 
open boundary conditions. Some additional details, which also apply to the 
simulations of the last section are given here. Although opening the boundary 
conditions might seem a rather drastic step, when larger lattices are used it is 
certainly practical and causes fewer difficulties than one might expect. The 
outer three layers were discarded in order to lessen the effects of the open 
boundary. The vast majority of the boundary-effect occurs here. In other 
words, for a 16 4 lattice quantities are only measured on the inner 10 4 and 
for finite-size scaling purposes the size is 10 4 . This can be thought of as a 
soft open boundary. All remaining boundary effects can be absorbed into 
the finite size scaling applied to different lattice sizes. In principle one does 
not need to exclude boundary layers, but that would require larger lattices 
for sensible results. One also loses the translational invariance, which, as 
mentioned above, prevents the calculation of the second moment correlation 
length. Otherwise analysis is relatively unchanged. 

Lattices of 16 4 , 20 4 and 24 4 were measured for 50,000 sweeps, after an 
initial equilibration of 10,000 sweeps, for (3 between 2.3 and 3.6 . The gauge 
overrelaxation was very slow at (5 = 2.3, taking roughly 4000 sweeps per 
Monte Carlo sweep. For (3 > 3.0, 100 or fewer overrelaxation sweeps were 
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Figure 8: Binder cumulant crossing for the 4-d SU(2) theory. 

needed. Overrelaxation was terminated when the per-link value of the gauge 
condition changed by less than 2 x 10 -8 . Test runs with a criterion ten 
times smaller showed no difference. Of course for many of these large (3 the 
correlation length is larger than the largest lattice size. Simulations designed 
to explore continuum physics must be run in a region where the correlation 
length is smaller than the lattice to avoid finite lattice and finite temperature 
effects. However, what is being sought here is a possible critical point, where 
the correlation length is infinite. So, just as in the study of magnetic theories, 
one must march right through the suspected critical region with finite lattices 
and interpret the data using finite-size scaling. If the true critical point is at 
(3 = oo then finite-size scaling should be able to see that too. In that case, 
the scaling in the entire range studied should look like the unbroken region 
with no apparent critical point. In particular, the Binder cumulant should 
everywhere be a decreasing function of lattice size. 

The Binder cumulant, U, shows an apparent crossing at quite a high 
cumulant value somewhere between (3 = 3.1 and 3.3 (fig. 8). The crossing 
around (3 = 3.2 is a bit messy, but for — 3.4 to 3.6 U is definitely an in- 
creasing function of lattice size indicating this is in the spontaneously broken 
region. Errors at these /3's are very small. At /3 = 3.4 U24 exceeds U20 by 4cr 
and Uie by ller. At (3 = 3.5 these are 5cr and 14cr respectively. Magnetization 
histograms are given in Fig. 9, showing a pattern typical of a higher order 
transition. For this 0(4) order parameter one must take into account the ge- 
ometrical factor (from solid angle) that biases the distribution toward larger 
magnitudes. In the unbroken phase, the distribution of magnetization mod- 
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Figure 9: Magnetization histograms at (3 = 2.3, 2.5, 2.9 and 3.5 on 24 4 
lattices. The latter three use the right-side vertical scale. 

uli, |m|, is expected to be a factor of |m| 3 times a Gaussian, exp(— |m| 2 /2<r^J. 
To more easily see the Gaussian behavior, the probability distribution P(|m|) 
is obtained by histogramming, and the quantity P(|m|)/|m| 3 is plotted. The 
value of |m| for each bin is not taken at the center, but at a value that would 
produce a flat histogram in an |m| 3 distribution, regardless of bin-size choice. 
This is 

H 3, n = 1 (Hl^?) (3) 
bm 4 (m 2 - mi) 

where m 2 and mi are the bin edges. This detail affects only the first couple of 
bins in the histograms. Magnetization and subtracted susceptibility curves 
are shown in Fig. 10 with the susceptibility showing a peak growing with 
lattice size and shifting toward weaker couplings. Scaling collapse fits were 
then performed for x an d < |m| > using U itself as the scaling variable. This 
can be done since U has no anomalous dimension, and has the advantage 
of being a single parameter fit (no fitting of /3 C is needed). From these fits, 
values of j/u = 2.951(3) and f3/v — 0.0250(8) were found. Since these are 
consistent with the hyperscaling relationship 7/1/ + 2/3/f = d = 3, it was 
decided to redo the fits with only a single degree of freedom by enforcing 
hyperscaling(Fig. 11), which yielded j/is = 2.950(2) with x 2 /d.f. = 1.7. 

Once 7/V has been determined, a correlation length can be defined from 
the bare susceptibility as above (Eqn. 2). Plots of should cross at an 
infinite lattice critical point. Such a plot is given in Fig. 12, which shows 
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Figure 11: Magnetization and susceptibility scaling collapse fits using U as 
the scaling variable. 
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Figure 12: Crossing plot for £l/L for 4-d SU(2). 

a rather clear crossing, also around (3 = 3.2. These crossings are verified 
by even larger confidence intervals than for U, with the 24 4 data exceeding 
the 20 4 by 18<r and the 16 4 by 32a at = 3.5. Finally a full four-quantity 
collapse fit for (3 C and v was performed utilizing U, x, < \m\ >, and £l/L 
(Fig. 13). This yields (3 C = 3.18(8) and v = 3.42(+0.32 - 0.23), with a 
X 2 /d-f = 3.2. Again, the 16 4 data were excluded from the fit. The slightly 
high x 2 /d-f- probably indicates the presence of some higher order corrections 
to scaling still operating at these lattice sizes. Errors in critical exponents 
and f3 c were determined by forcing the given quantity higher or lower until the 
X 2 /d.f. doubled, with the other parameters free to change. This somewhat 
more conservative method than "adding unity" to x 2 jd.f. guards against 
any possible underestimate of statistical errors, which were from plateaus of 
binned fluctuations. 



5 Fits to gauge invariant quantities 

When one thinks of confinement, the quantity that immediately comes to 
mind is the string tension. In particular, the inverse square root of the string 
tension defines a correlation length. If the phase transition seen above is 
to be consistent, the correlation length in the confining phase should match 
that determined from the string tension. If the string tension scales with 
the same exponent given above, that would help corroborate the existence 
of the phase transition, since the string tension is determined from gauge 
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Figure 13: SU(2) scaling collapse plots for U, x, < 1^1 >> an d £l/L. 
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invariant Wilson loops and does not involve any sort of gauge fixing. One 
can also attempt fit the "singular part" of the specific heat in this region to 
the expected scaling behavior. Finally, the positions of the deconfinement 
transition observed on lattices with one short dimension, usually interpreted 
as a finite temperature transition can be analyzed as finite-lattice /3 -shifts of 
the infinite lattice transition seen above, which also depend on the exponent 
v. It will be shown that all three of these gauge invariant quantities scale 
consistently with the above transition seen at (3 C = 3.18 with v = 3.4. 

However, first the three-dimensional scaling of the hyperlayers needs to be 
reconciled with the four-dimensionality of the original system. For instance 
what value of d should be used in the hyperscaling relationship for the specific 
heat exponent a = 2 — dpi The key is to consider the relationship between 
the correlation length being measured from the spin correlations in the 3- 
d hyperlayers and that of the four-dimensional theory (being measured by 
the string tension for instance). The magnetization correlation function can 
be considered similar to the lxN Wilson loop or, more accurately, a unit 
length partial Polyakov loop(PPL) or "gaugeonjTj." In the Coulomb gauge a 
segment of links of any length, aT in the fourth direction is an observable. 
Here T is the time extent in units of the lattice spacing. The correlation 
function between two PPL's a distance Ra apart is expected to behave as 

exp(-aRTa 2 ) (4) 

where a becomes the usual asymptotic string tension for large T[T]. For 
T = 1, a will differ somewhat but the dimensionality will still be the same. 
The a for T = 1, sometimes referred to as the Coulomb string tension, has 
been observed to approximately scale with the asymptotic string tension 
[20] . It is therefore an inverse length-squared object. However, this same 
object was treated as the inverse correlation length itself in the 3-d analysis, 
because there the 4-d links were treated as "spins" with no inherent length 
associated with them, whereas in the 4-d theory they are associated with the 
lattice spacing. Therefore the correlation length determined from the spin 
analysis is really an area in the 4-d theory. Thus the correct 4-dimensional 
exponent for a true length, 1 / y/a, should be half of the 3-d measured value 
of 3.42, i.e. i/ 4 = 1.71. 

Neither the string tension nor the specific heat show much finite-size vari- 
ation in the /3 region in which they have been measured (the region where the 
correlation length is less than L). Therefore, they can be expected to rather 
closely follow infinite-lattice scaling laws. Fig. 14 shows a plot of V ad 2 vs. 
j3 with a one-parameter fit to the form c(l//3 — l/(3 c ) Ui . String tensions were 
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Figure 14: String tension fit to finite-order scaling law. 



taken from refs 0|2T]. A somewhat reasonable, though not entirely satisfac- 
tory, fit is obtained for the values fi c = 3.18, = 1.71 obtained above from 
the 3-d magnetization. The fit particularly differs with the (3 = 2.85 point(6 
standard deviations off). Without this point the x 2 /d-f- — 3-5. If v is allowed 
to vary, only a slightly better fit is obtained for z/ 4 = 1.67 (x 2 /d.f. = 4.08 
- fit improvement does not compensate for increasing d.f. by 1). There- 
fore the string tension appears to be selecting essentially the same value for 
the critical exponent as obtained from the Coulomb gauge magnetization. 
Since string tensions are measured with a variety of methods and on differ- 
ent lattice sizes, some systematic differences between the datapoints may be 
present, which could explain some of the difficulties with the fit. Similar fits 
can be made for a fairly broad range of f3 c , so the string tension is not good 
at pinning that down further. 

The specific heat exponent is expected from hyperscaling to be a = 
2 — 4z/4 = —4.84. The high negative value indicates a very soft singular- 
ity. Five derivatives of the specific heat would have to be taken before an 
infinite singularity would be seen. Such a weak singularity would be virtually 
impossible to detect directly using numerical methods and certainly explains 
how this transition could have been missed before. Such high values of v are 
unusual in spin models, but more common in spin glasses, because the lower 
critical dimension (led) of many spin glasses is thought to lie around 2.5 [22J, 
and v — > oo at the led. Since it is 3-d hyperlayers that are being studied 
here, and a spin-glass to ferromagnet transition is being hypothesized, the 
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Figure 15: Specific heat finite-order scaling fit. Lower curve is non-singular 
part of fit. 



closeness of the dimensionality to the led could explain the unusually high 
value of v. For a finite singularity, scaling can still be checked. Since the 
finite size dependence of the specific heat is known to be very small, one 
can fit to the expected infinite lattice scaling function. The following fitting 
function was used: 

C {l/P - l/f3 c )^-V + d + C 2 /f3 + C 3 /f3 2 + Q//? 3 . (5) 

The last four terms are needed to fit the non-singular part. C\ is set to 
the perturbative value of 0.75, but C2 through C4 are allowed to differ from 
the perturbative values to account for higher-order neglected terms. Fig. 15 
shows that data in the range 2.44 < (3 < 3.2 fits well to this form with the 
previously determined v± = 1.71 (x 2 /d.f. = 0.58). Data were from plaquette 
fluctuations on 24 4 lattice for runs of 500,000 sweeps or more. Again if z/4 is 
allowed to vary, the best fit is found at a nearby value of 1.63 (x 2 /d.f. = 0.57). 
Without the 1//3 3 term a value of 1.53 is preferred (x 2 /d.f. = 0.62) The 
specific heat is not very sensitive to the exact value of (3 C . 

Finally one can look at the scaling of (3 C on lattices with one short dimen- 
sion - what has usually been interpreted as a finite temperature deconfine- 
ment transition. If there is an infinite-lattice critical point at (3 = 3.2, then 
the phase transition on the distorted lattice would not be a finite-temperature 
transition, since the continuum limit is already deconfined at zero tempera- 
ture. Rather it is simply the same transition /3-shifted by the finite size shift 
effect. The shift in the critical point for a finite lattice is expected to be[i~8] 

(& - (3 CT ) oc L- 1 '^ (6) 
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Figure 16: Scaling fit for finite-temperature transition points according to 
the finite-order zero-temperature phase transition hypothesis(straight lines). 
Diamonds are for the Wilson action and squares are a Z2 monopole sup- 
pressed action. Curves are the two-loop renormalization group prediction 
of the no phase transition hypothesis, made to pass through the highest L T 
points. 

for L T « L, where L T is the temporal lattice size. For large v this can be 
substantial and dies only slowly with lattice size. If one plots (3 CT vs. L^ 1 ^ 4 , 
with z/4 = 1.71 (Fig. 16), a reasonable linear fit is seen to Wilson-action data 
from Ref. [6], with an intercept giving (3 C = 3.11(4), close to the previously 
identified value of 3.18(8). Also shown is the two- loop weak-coupling renor- 
malization group prediction for the no phase transition hypothesis, which 
does not fit the data as well. In addition the same quantity but for an al- 
ternative action which suppresses Z 2 monopoles, as studied by Gavai[2"3"]. is 
plotted. These data also fit well to a linear fit with an intercept /3 C = 3.13(4), 
in agreement with the Wilson value, and do not fit at all well with the per- 
turbative renormalization group. The monopole suppressed action has the 
same weak coupling perturbation series as the Wilson action. For weak cou- 
plings, results should merge with the Wilson action ones, so the (3 C values 
would be expected to be close, though not necessarily identical. The system 
with L — > 00 with L T finite is, of course, a true 3-d system and thus will have 
different critical exponents from the 4-d theory. As L T — > L, a dimensional 
crossover is to be expected with critical exponents morphing to their 4-d 
values [21]. 
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Therefore, the gauge invariant quantities of string tension, specific heat, 
and finite-lattice shift all scale reasonably in accordance with the zero tem- 
perature phase transition seen by the Coulomb gauge magnetization order 
parameter. This lends further credence to the existence of this phase transi- 
tion - in particular that gauge fixing is not necessary to see it. It is merely 
a very convenient and probably necessary step to take in order to define a 
local order parameter. 

The idea that the pure gauge SU(2) lattice gauge theory might have a 
lattice-artifact driven phase transition in analogy to U(l) has been consid- 
ered before [251 126]. Perhaps the best previous evidence to date that such a 
transition must exist involves another expansion of the coupling space, the 
fundamental- adjoint plane. For a sufficiently large adjoint coupling a first 
order bulk transition appears in the SU(2) theory. The scaling of latent heat 
with the size of the hysteresis region was shown in Ref. [27J to be inconsis- 
tent with the endpoint of the first order line being an ordinary critical point. 
Rather it is only consistent with being a tri-critical point associated with 
an exactly broken symmetry, for which the transition cannot just end, but 
must continue as a higher-order one, again until it hits an edge of the phase 
diagram. Slopes of lines of constant physics predict a crossing of the Wilson 
axis. Consistent with this hypothesis, when "finite temperature" transition 
points are followed into the fundamental-adjoint plane they appear to head 
toward this tricritical point [SB]. Exact joining is difficult to establish, how- 
ever, and has been disputed|2§], but if the endpoint of the first-order line is 
truly a tricritical point then joining must take place. 

6 Conclusion 

It has often proven useful to connect theories by using an extended coupling 
space, especially when the behavior of one theory is better known than the 
other. This paper extends an old suggestion by Darhuus and Frohlich relat- 
ing the 4-d SU(2) lattice gauge theory to the 3-d 0(4) Heisenberg modelfTTj. 
by allowing the couplings of horizontal and vertical plaquettes to differ. The 
enlarged model becomes the 3-d 0(4) Heisenberg model in the — > oo limit 
and the 4-d SU(2) lattice gauge theory in the (3h = fiv limit. A Monte Carlo 
simulation at 0h = 20 shows a phase transition at /3y c = 1.01. In order to 
follow the same ferromagnetic order parameter as in the spin model, the lat- 
tice gauge theory must be transformed into Coulomb gauge, in which the 4th 
direction pointing links act as 0(4) spins exhibiting a 3-d global symmetry, 
the remnant symmetry of the gauge fixing. Although Coulomb gauge fixing 
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has been plagued with large systematic and random errors introduced by 
the gauge fixing algorithm, it is found that using open boundary conditions 
solves these problems - with errors hundreds of times less than when using 
periodic boundary conditions, nearly eliminating the lattice Gribov problem. 
The local order parameter allows standard methods for finding phase tran- 
sitions such as Binder cumulant crossings and finite-size scaling collapse fits. 
Not surprisingly, because of the high value of (3h, the flu = 20 transition 
is very similar to the Heisenberg transition at (3h = oo, (3y = 0.936, with 
similar critical exponents. However, the fact that the transition still exists 
for non-infinite (3h has far-reaching implications. The ferromagnetic transi- 
tion breaks an exact symmetry, the remnant gauge symmetry left over after 
Coulomb gauge fixing. Therefore, it must divide the coupling plane into two 
distinct regions - symmetry broken and unbroken. The line of phase tran- 
sitions must persist until it hits another edge. The Heisenberg model itself 
has only one critical point, so the line of phase transitions must cross the 
line 0v = & at a finite value (somewhere between 1.01 and 20 - as (3h de- 
creases, fiy c increases). Indeed a high-order transition is seen around (3 = 3.2, 
showing Binder cumulant and crossings, and good collapse of various 

graphs to scaling functions about this point. The transition also predicts 
scaling laws consistent with the behavior of the specific heat, string-tension, 
and even finite-temperature transition points (reinterpreted as finite-lattice 
shifts) in the crossover region. 

Much of the lattice gauge theory program has been based on the assump- 
tion that the non-abelian gauge theories SU(2) and SU(3) have no phase 
transition and the continuum limit is confining. Of course the SU(3) case 
needs to be checked in detail, but if SU(2) has a zero-temperature phase 
transition it is likely SU(3) has one as well. A non-confining and symmetry- 
broken continuum limit would require a new mechanism of quark confine- 
ment, since it would no longer be a property of the gluon sector acting alone. 
One possibility is that the light-quark chiral condensate is itself the source 
of confinement, a hypothesis which has been advanced by Gribov [30J and 
others [3Tj. A strong enough force, whether confining or not can cause quark- 
antiquark pairs to condense in the vacuum, breaking chiral symmetry. If this 
chiral vacuum expels strong color fields, of which there is some evidence[32j, 
then the chiral vacuum itself could cause a bag-like pressure around strong 
color sources, with the lowest energy configuration of vacuum + hadron being 
one in which the color fields are confined to a small volume. In most cases 
unbreaking of chiral symmetry from thermal effects and deconfinement are 
coincident, which supports such a linkage. Thus it is possible that simulations 
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with light quarks are exploring correct continuum physics (whereas the pure- 
glue simulations in the confining phase are exploring a strong-coupling phase 
not connected to the continuum limit). To be sure, alternative gluon actions 
need to be explored which eliminate, if possible whatever lattice artifacts 
are causing the pure-glue phase transition, in analogy with the monopoles 
of the U(l) theory. Preliminary data suggest that eliminating both negative 
plaquettes (Z2 vortices and monopoles) and a different kind of monopole de- 
fined from the non-abelian Bianchi identity, which I have called SO(3)-Z2 
monopoles [33], may force the SU(2) theory to remain in the ferromagnetic 
phase for all couplings. 

The symmetry breaking which occurs in the weak coupling phase is not 
unexpected. The U(l) theory in the Coulomb gauge almost certainly is 
also magnetized in the same way, becoming paramagnetic as the coupling 
becomes stronger at the monopole-induced phase transition to a confining 
theory, but being magnetized in the weak-coupling region down to the weak- 
coupling continuum limit. In this gauge, at least, the photon is a spin- wave 
like Goldstone boson associated with the breaking of the remnant symmetry 
left over after Coulomb gauge fixing, as opposed to the original bare gauge 
field the theory started out with. This picture has been previously explored 
in the continuum|34j. The Gupta-Bleuer quantization method may be an 
alternative way to see how the elementary fields build into a Goldstone boson, 
and how a linked gauge and Lorentz symmetry is "reborn." A similar picture 
may hold for the pure-gauge SU(2) and SU(3) gauge theories, with the gauge 
fields remaining massless and the running coupling having an infrared fixed 
point. However, when light quarks are added, this vacuum could become 
unstable and confinement may emerge as a byproduct of chiral symmetry 
breaking, following Gribov's or a similar mechanism. 
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